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ENERGY  TRANSFER  FROM  HIGH  POWER  PULSELINES 
TO  THE  NEXT  GENERATION  OF  PRS  LOADS 
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I.  lutroductiou 
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'  The  design  of  more  complex  loads  for  Plasma  Radiation  Sources  and  the  quest 
for  ever  higher  simulator  radiative  yijeJds  and  power  levels  require  heretofore  unex¬ 
plored  energy  densities  in  the  load  region,  perhaps  as  high  as  1.50  Quite 

apart  from  detailed  design  questions,  there  is  ^ome  concern  that  these  higher  power 
and  energy  densities  (delivered  by  ^sub-Ohm^  pulselines)  cannot  couple  well  to  the 
PRS  configurations  used  today^^^ir  to  even  higher  inductance  designs  that  might  he 
used  in  the  future. 

In  this  report  we  will  examine  the  energy  transfer  to  two  PRS  loads  from  generic 
pulselines  of  the  Saturn  class  operating  at  the  8  — 15  MA  level.  Saturn  is  nearly 
a  1  Ox  machine  in  terms  of  power  flows  and  has  recently  been  outfitted  with  PRS 
front  end  hardware;  it  is  modeled  in  a  .subsequent  section  (III).  A  basic  slug  PRS 
model  and  an  appropriate  X-pinch  model,  derived  from  2-D  MHD  simulations,  are 
used. 


The  pulselines  are  modeled  with  a  transmission  line  code,  described  in  the 
next  section  (II).  Apart  from  generalization  to  spacetime  dependent  line  parame¬ 
ters  (impedance,  propagation  speed,  and  damping),  the  solution  technique  is  quite 
standard,  robust,  and  accurate  in  its  energy  transfer  characteristics.  Even  the  .sim¬ 
plest  transmission  line  models  of  machines  like  Saturn  provide  the  expected  result 
-  PRS  loads  will  accommodate  the  line  impedance  and  draw  energy  effectively.  A 
more  delicate  question  is  that  of  loss  mechanisms.  This  is  left  to  future  work. 

In  Section  IV  the  coupling  of  pulsers  like  Saturn  or  the  Falcon  design  (proposed 
by  PI  as  a  lOx  device)  to  an  X-pinch  load  is  discussed.  The  particular  load  chosen 
is  first  modeled  with  the  2-D  MHD  code  PRISM  to  exhibit  the  behavior  we  wish  to 
exploit.  The  inductance  history  can  then  be  u.seil  as  a  rough  guide  to  the  behavior 
Manmcript  approved  November  29,  1988. 
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of  such  loads  in  Saturn  or  Falcon  class  machines. 

Section  V  is  devoted  to  a  discussion  of  the  results,  conclusions  and  the  impli¬ 
cations  for  future  efforts. 

11.  l^ausnussiou  Line  Model 

Existing  transmission  line  models  come  in  several  variations.  Within  the 
MAGIC  ^  program  exist  options  for  coupling  transmission  lines  to  a  variety  of  2-D 
simulation  regions  as  well  as  the  use  of  parallel  and  series  junctions  among  lines. 
Within  the  BERTHA  *  package  exists  a  broad  spectrum  of  line  elements  on  which 
waves  are  propagated  without  distortion  by  summation  of  forward  and  backward 
disturbances.  In  any  study  of  energy  coupling  to  the  time  dependent  PRS  one  must 
join  the  transmission  line  to  varying  load  parameters  Z(t),  L(R(t)),  allow  for  space 
time  varying  line  parameters,  and  also  drive  the  system  with  a  fixed  initial  energy. 
No  existing  tratismission  line  code  offered  ready  access  to  all  these  capabilities,  so  a 
simple  implementation  of  the  1-D  telegrapher's  equations  has  been  written  to  treat 
such  problems. 

Resolving  the  continuous  transmission  line  into  discrete  units  of  inductance  per 
unit  length  (L')  in  series  with  the  load  and  capacitance  per  unit  length  (C')  in  par¬ 
allel  with  the  load,  the  telegrapher's  equations  can  be  written  in  the  dimensionless 
forms 


nAly,  ^  ny.  _  ("+1/2/.  _  "+l/2/^_j)^ 


(T)_  f  ("Vy+,  -  "Vj)  , 


where  the  matrix  operator 


introduces  an  adjustable  damping  parameter 


1  kN 

-  oc  - — 
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mp 


which  controls  the  dispersion  of  the  method  at  high  wavenumber.  By  damping  out 
the  large  wave  number  components,  the  tendency  of  these  difference  equations  to 
propagate  higher  wavenumbers  more  slowly  is  mitigated  and  low  wavenumber  pulse 
shapes  are  then  transmitted  with  low  harmonic  distortion.  Here  the  parameter 
a  =  cht/hx  is  the  Courant  number,  and  Zj  can  be  spatially  and  temporally  variable. 

At  the  boundaries  the  damping  is  set  to  zero,  the  combination  of  exterior  circuit 
elements  and  the  local  currents  on  the  transmission  line  then  fixes  the  time  advanced 
currents  and  voltages  through: 


n+l/2  j 


> 


tl+1 


/Vy  =  rVj+i  -  "K,)  , 


”V 

'  o 


ZotVM)  +  KVi 

ZqT  +  ho 


rZj  ^IjZL(t)  +  LjRrVj 
tZj  +  L{R) 


where  the  time  derivative  of  current  has  been  eliminated  in  favor  of  the  appropriate 
voltage  differences  and  r  is  the  delay  interval  of  a  single  element,  viz.  r  =  h^/c. 
The  source  voltage  V„{t)  can  he  derived  through  a  variety  of  methods,  e.g.,  a  fixed 
waveform,  or  a  computed  capacitive  discharge.  The  load  is  characterized  by  a  time 
dependent  inductance  L{R,t)  and  resistance  Z[,(t).  Here  R  denotes  the  load  radius, 
and  R  must  be  included  in  the  load  resistance. 

For  a  PRS  load  a  wide  variety^  of  simple  models  has  been  used  in  this  context. 

Each  choice  made  in  configuring  the  load  model  has  consequences  for  the  calcu- 
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Jation  as  a  whole  which  are  substantially  impossible  to  generalize  and  which  pro¬ 
duce  a  variety  of  transmission  line  responses.  In  order  to  set  the  proper  damping 
parameter  and  thus  to  calibrate  the  basic  method,  a  simpler  problem  is  appropriate. 

If  the  input  voltage  is  specified  as  a  Gaussian  pulse  with  a  variable  width  rela¬ 
tive  to  the  transmission  line  time  and  the  load  is  a  fixed  impedance  equal  to  that  of 
the  transmission  line,  then  the  ideal  response  is  a  distortion  free  propagation  of  the 
Gaussian  pulse.  Given  an  upper  limit  to  the  wavenumber  which  can  he  resolved  by 
the  difference  equations,  the  damping  added  must  attenuate  the  higher  wavenumber 
components  enough  to  prevent  the  ’’shedding”  of  aphysically  slow  ’’wake”  distur¬ 
bances,  but  not  so  large  as  to  interfere  with  an  accurate  portrayal  of  the  energy 
transport  in  the  line. 

The  use  of  a  7  value  of  about  0.125,  corresponding  to  a  damping  wavenumber 
at  about  4  times  the  Nyquist  wavenumber,  offers  a  reasonable  compromise  in  this 
tradeoff.  With  an  overall  energy  conservation  error  of  5.1 5  *10”^,  Figure  1  shows  the 
result  of  propagating  a  50  ns  pulse  over  a  22.5  ns  length  oflfSfl  line  into  a  matched 
load.  The  Gaussian  envelope  at  the  output  is  essentially  self-similar  to  that  at  the 
input,  being  scaled  down  slightly  due  to  the  attenuation  in  the  method.  There  are 
no  refiections  at  the  interface  for  the  matchload  at  any  part  of  the  wavenumber 
spectrum. 

Figure  1 
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III.  Aoplicatiou  to  Saiuru 


A  simple  model  of  the  Saturn  machine  at  SNL  can  be  formed  from  three  ele¬ 
ments:  a  pulse  source,  a  long  transmission  line  section  equivalent  to  the  four  parallel 
water  lines  which  bring  the  power  down  to  a  radius  of  1  m,  and  a  short  transmission 
line  delivering  the  power  to  the  load  region  which,  in  the  PRS  mode,  offers  an  initial 
inductance  of  5— *10  iiH. 

In  practice  the  transition  between  the  two  lines  is  a  post-hole  convolute,  while 
the  shorter  line  offers  a  fixed  impedance  for  only  part  of  the  path  to  the  load,  after 
which  it  becomes  a  simple  radial  line,  viz.  Z  oc  1/i?.  The  details  of  the  Marx  bank 
will  be  ignored  here  -  let  it  be  abstracted  to  a  single  fast  pulse  (4.4  MV  and  50  ns  ^ 
width)  in  series  with  a  small  resistance  (0.12  ft).  Likewise,  the  convolute  details 
will  be  ignored  as  well,  it  is  abstracted  to  a  discrete  jump  in  impedance  to  3.3  ft, 
lip  from  the  1/4  ft  of  the  parallel  water  lines. 

In  terms  of  detailed  experimental  support  such  minutiae  as  these  are  important, 
but  in  this  examination  of  gross  power  coupling  they  will  not  matter  and  they  are 
modeled  *  elsewhere.  The  primary  quest  here  is  the  energy  coupling  characteristics 
of  Saturn  class  machines  to  various  PRS  configurations.  The  simplest  adequate 
model  can  therefore  be  abstracted  to  the  following  transmission  line  problem. 

Figure  2 
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To  drive  lOx  machine  currents  to  a  PRS  load  in  %  lOOna  requires  Voc  ^  AMV  on 
the  driver  pulse,  and  a  forward  energy  Auence  of  «  IMJ  into  the  waterline. 

As  a  further  benchmark  test  of  the  transmission  line  code,  a  similar  problem 
was  run  on  the  BERTHA  package^^  with  an  1/8  fl  line  impedance  and  a  Axed 
load  impedance  of  3 /S  Q.  Both  calculations  agree  to  a  few  percent  in  the  timing 
of  various  reAections,  the  timing  of  current  reversals  and  in  the  peak  values  of  the 
voltage  and  current  waveforms  (at  the  input  to  the  line  and  at  the  load). 


Slug  PRS  Model 


The  early  phase  of  implosion  presumes  a  snowplow  like  compression  of  the  full 
load  mass.  Even  though  such  an  idealization  is  not  likely  to  he  completely  accurate, 
for  several  reasons,  the  most  serious  consequence  is  a  shift  of  the  optimum  mass 
point  for  radiation  yield.  The  most  important  variable  in  the  calculation  of  the  run 
down  is  the  motional  impedance.  I  have  shown  elsewhere^  ^  that  a  variety  of  details 
in  the  current  penetration  and  radiation  physics  still  allows  a  rough  adherence  to 
the  familiar  snowplow  velocity  scaling  with  electric  held  -  either  for  a  slug  or  for 
a  fully  resolved  1-D  load  model.  As  a  consequence  the  slug  model  will  capture  the 
proper  trends  in  generator/load  coupling  even  if  the  point  by  point  comparisons 
may  differ  somewhat  from  the  actual  experiment. 

Employing  a  series  resistance  and  inductance  for  the  slug  load,  both  time  vary¬ 
ing  as  determined  by  the  load  radius,  the  governing  equations  for  the  run  down 
phase  become: 


JL. 

pc^R  ' 


T  (  rj\  I  /  > 

L(H)  =  —  ln(— , 
c*  n 


Zl 


Zo- 


UR 

c2  R  ' 
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with  Zo  fixed,  d  the  AK  gap,  Rwaii  the  radius  of  the  return  current  path  in  the 
load  region,  p  the  load  mass  per  unit  length,  and  Zg  the  average  impedance  of  the 
plasma,  viz.  Zo  =  The  choice  rj  is  reasonably  taken  to  produce  a  Zg  in  the 

10m fl  range  at  early  times.  Coupling  the  relations  to  the  transmission  line  solver 
is  accomplished  by  calculating  R,  R,  Z^  and  L(R)  at  each  timestep.  The  energy 
absorbed  in  the  load  region  is  then  J*  dt\  Vi^{ti  as  inferred  from  the  line 

solution. 

As  the  run-down  terminates  the  plasma  load  should  ideally  undergo  a  non- 
adiabatic  stagnation  and  convert  the  input  energy  to  radiation.  The  final  state 
of  the  load  plasma  is  almost  never  very  simple,  but  a  rough  average  picture  is  to 
model  it  as  a  Rennet  equilibrium  pinch.  The  Rennet  current  figures  as  a  central 
notion  in  oscillatory  MHD  solutions  as  well,  so  an  order  of  magnitude  sort  of  model 
can  set  the  final  impedance  at  a  fixed  final  radius,  Rmin-  If  the  pinch  temperature 
is  eliminated  in  favor  of  the  current,  oii"  finds  that  the  Rennet  profile  offers  an 
impedance 


Zl 


1.109  - 10®( 


d 

02  rs 

^min‘kA 


)|n| 


where  is  in  units  of  kA,  is  chosen  initially,  and  the  motional  impedance 

is  assumed  to  vanish. 


Energy  Coupling  and  Stagnation  Radius 

The  radial  structure  of  a  PRS  load  will  generally  have  a  substantial  effect  on 
the  radius  of  the  current  path  at  stagnation.  Efficient  radiative  loss  can  deepen  the 
implosion  by  perhaps  a  factor  of  10:  turbulent  enhancement  of  plasma  resistivity 
can  spread  the  current  channel  and  increase  the  final  radius.  While  the  final  radius 
is  therefore  not  a  "free"  parameter  is  any  sense,  it  is  a  very  relevant  measure  of  the 
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implosion  quality.  Moreover,  the  logarithmic  singularity  in  the  load  inductance  as 
Rmin  —*  0  makes  this  parameter  appear  to  be  important  for  energy  transfer.  While 
the  values  of  Rmin  examined  here  are  a  bit  smaller  than  those  observed  in  practice, 
the  machine  behavior  expected  if  very  small  radii  are  achieved  is  a  worthwhile  target 
in  a  study  of  this  kind. 

On  the  transmission  line,  the  load  implosion  and  stagnation  appears  as  a  back¬ 
ward  propagating  pulse  -  lowering  the  current  and  increasing  the  voltage.  As  the 
load  samples  larger  inductance  values  due  to  smaller  Rmin  ^his  ’’downward'’  current 
pulse  turns  into  a  full  current  reversal.  The  larger  Rmin  values  thus  imply  a  mild 
power  reflection  back  up  the  line,  while  the  smaller  values  imply  a  rather  effective 
rejection  of  power  by  the  load. 


Figure  3(a,b) 


Vi  at:  92.5  |ns| 

I  I . I,.  -I 

^  . 


<  I  ■  s  •  <  • 

02  04  00  00  10 


Figure  4(a,b) 


-0.5 


VJ  at:  too  M 

I  I  -  * 


U\ 


04  06  08  10 


0  25.0  SO.O  7S.0  100.0  125  0n» 

LOAOKE 
502224  2501 


8 


In  Figure  3  are  two  snapshots  of  the  reffection  in  a  Saturn  calculation,  one  at 
Rmin  =  0.01cm,  the  other  at  Rmin  =  O.OOlcm.  In  Figure  4  are  shown  the  drive 
and  load  I  and  V  waveforms,  as  well  as  the  energy  into  the  line  and  into  the  load. 
The  load  radius  is  also  plotted.  All  such  runs  as  these  were  done  with  a  load  mass 
of  100  pg/cm  and  produced  a  peak  load  current  ^  8MA.  The  mass  is  such  that 
implosion  occurs  at  ss  lOOn^  and  the  rejected  power  is  therefore  unable  to  go  up 
and  back  before  the  run  ends  at  125  ns. 

Under  the  assumptions  of  this  PRS  model,  when  the  load  reverses  the  current 
the  energy  leaves  the  load  region  because  all  the  energy  is  in  the  fields  around  it. 
The  real  behavior  would  likely  be  more  dissipative  and  not  allow  as  much  energy 
release  to  the  line.  It  is  therefore  reasonable  to  compare  the  energy  coupled  to  the 
load  region  to  that  coupled  to  a  1/4  ft  matched  load  in  the  constant  125  ns  interval. 
The  results  are  summarized  in  the  following  Figure  5,  where  the  efficiencies  of 
overall  energy  transfer  -  =  Eioaj/ Ematch  -  transfer  to  load  kinetic  energy 

-  rjT  =  TpRs! Emaieh  afc  plotted. 

Figure  5 
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IV.  X-pinciies  in  Saturn  and  fVilcon 

The  X-pinch  Joad  Jjas  heretofore  been  configured^ as  two  or  more  crossed 
wires  in  a  diode  gap.  The  central  motivation  being  to  focus  more  of  the  machine 
power  on  a  smaller  volume  of  plasma  in  hopes  of  exciting  more  K-shell  radiation. 
A  different  quest  for  the  X-pinch  can  also  be  contrived  -  the  production  of  a  "mini¬ 
diode"  plasma  state  which  will  make  effective  use  of  nonthermal  electron  populations 
in  exciting  K-shell  radiation. 

In  order  to  produce  such  a  "mini-diode"  one  can  exploit  the  sausage  instability 
through  judicious  mass  loading.  Using  the  simple  scaling  law  for  the  implosion  time 
pointed  out  by  Apruzese  ^ ,  a  constraint  of  fixed  time  from  different  radii  imposes  a 
particular  mass  per  unit  length  on  any  extended  load  configuration  - 

=  1.1972. 10’(^)-  Ipg/cmj, 

such  that,  for  a  fixed  run  down  time  Tg  to  a  peak  current  Ip,  a  figure  of  revolution 
RoiZ)  must  be  mass  loaded  as  above  in  order  to  implode  synchronously  at  all  axial 
locations.  In  Figure  0  are  shown  the  density  contours  of  such  a  load  exactly  as  it 
was  represented  in  the  2-D  MUD  code  PRISM.  The  conical  portion  helps  to  keep 
the  inductance  lower  at  early  times,  while  the  cylindrical  segment  near  the  midplane 
is  allowed  a  greater  or  lesser  rundown  length  to  provide  a  store  of  kinetic  energy 
upon  .stagnation. 

Figure  0 
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This  Joad  configuration  was  imploded  through  the  use  of  a  current  ramp  which 
rose  linearly  to  3  MA  in  100  ns.  While  there  is  certainly  some  detailed  structure 
produced  during  the  rundown,  the  dominant  kinematics  is  precisely  as  inferred  from 
the  0-D  scaling  law.  By  90  ns  (c.f.  Figure  7),  the  conical  portion  has  run  down  a 
bit  faster  and  two  density  peaks  have  appeared  about  the  midplane.  Some  material 


Figure  7,8,9 
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has  been  pushed  ahead  near  the  midplane  due  to  magnetic  field  penetration  and 
the  resulting  spallation  when  the  shock  emerged  into  the  interior. 

Upon  stagnation  at  about  100  ns  the  load  dynamics  does  reflect  the  structural 
details  impressed  during  the  rundown.  The  sausage  instability  goes  rapidly  to  a 
saturated  state  -  the  current  density  is  shown  to  leave  the  necked  down  region  - 
and  the  originally  cylindrical  segment  is  bifurcated,  c.f.  Figure  8.  The  split  state 
persists  for  some  10  ns  more  through  the  end  of  the  simulation,  Figure  9. 

Here  ’’stagnation*'  connotes  the  arrival  of  material  on  axis,  the  history  of  the 
density  there  is  shown  in  Figure  10.  The  overall  load  design  has  thus  produced  its 
intended  result  -  there  is  a  density  cavity  of  nearly  two  orders  of  magnitude  near 
the  midplane  over  which  the  diode  voltage  must  be  dropped! 

Figure  10 

DENSITY 
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We  leave  to  future  work  the  assessment  of  runaway  electron  populations  in 
such  a  load.  For  now,  in  the  context  of  a  lOx  machine,  it  is  reasonable  to  examine 
the  energy  delivery  to  such  a  load  when  it  is  characterized  by  a  time  dependent 
inductance  and  the  implied  motional  impedance.  When  the  evolution  in  2-D  shown 
above  is  reduced  to  an  inductance  history  one  hnds  the  behavior  shown  in  Figure 
11.  There  is  shown  the  inductance  added  by  the  imploding  X-pinch  and  a  similar 
(in  mass,  radius  and  current)  slug  implosion. 

For  the  model  to  work  with  the  transmission  line  code  this  function  must  be 
mapped  into  the  familiar  triad  (  L(R),Zi,{R,  R),  and  R(t)  ].  The  radius  we  use  is  a 
current  weighted  one,  viz.  derived  from  the  volume  exterior  to  the  current  channel. 

Figure  11 
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Saturn  X-piuch 

When  the  X-pinch  model  is  used  in  place  of  the  slug  in  the  same  machine 
siniuJation  of  Saturn  detailed  above,  the  energy  transfer  is  somewhat  lower  due 
to  the  lack  of  enough  motional  impedance  to  draw  energy  from  the  line.  The 
overall  waveforms  are  shown  in  Figure  12,  and  the  inductive  ’’notch”  is  much  less 
pronounced.  The  inferred  efRciency  r]i  is  0.45,  and  the  peak  load  current  is  8.3M A. 

Figure  12 
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Since  the  Falcon  design  is  concerned  with  an  inductive  storage  device,  the 
circuit  and  transmission  line  model  changes  slightly.  In  a  Falcon  lOx  machine  eight 
inductive  energy  stores  (lES)  are  switched  together  into  a  single  MITL  at  about  a 
meter  radius  which  then  must  converge  the  power  to  a  PRS  load.  The  lES  modules 
would  have  a  quarter  period  of  %  1.27ps,  each  being  charged  to  perhaps  720  kV  and 
using  a  storage  inductor  of  about  150  nil  in  vacuum  prior  to  the  opening  switch.  A 
reasonably  generous  estimate  of  the  inductance  downstream  of  the  opening  switch 
(before  the  PRS  load  implodes)  is  pe-f  aps  18  nff.  Thus,  working  from  a  meter 
radius  or  3  ns,  the  final  transmission  me  impedance  is  6  H. 

When  all  eight  modules  are  cm  lected  and  switched  out  near  peak  current  one 
has  some  30  MA  available  *  j  drive  the  final  line  and  PRS.  Not  all  of  this  current 
is  expected  in  the  imploding  load,  most  of  the  energy  will  be  dissipated  in  the 
opening  switch  or  left  in  the  inductance  of  the  store  and  final  transmission  line. 
However,  even  with  all  these  sinks  the  final  energy  to  an  imploding  load  is  easily 
competitive  with  Saturn.  If,  for  example  the  Falcon  opening  switch  were  to  achieve 
an  R  K  O.OSn/ns,  the  energy  delivery  efficiency  to  a  slug  load  (qi, ),  relative  to  the 
original  stored  energy,  is  ss  0.46  -  in  the  same  range  that  Saturn  achieves  into  a 
slug  load.  This  efficiency  is  computed  over  the  same  time  interval,  [0, 125]ns,  and 
decays  with  slower  opening  rates.  If  the  X-pinch  load  is  examined  at  the  same 
opening  rate,  the  efficiency  decays  to  %  0.36.  The  peak  current  expected  in  the 
Falcon  X-pinch  would  he  «  17.QMA,  again  comparable  to  Saturn. 

V.  Couciusious 

All  the  experience  to  date  with  this  transmis.sion  line  and  load  model  indicates 
that  for  all  of  the  generic  machine  designs  presently  contemplated  for  the  lOx  ma- 
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chine,  for  aii  the  likely  domains  of  load  stagnation  radius,  one  can  state  two  results 
with  some  confidence. 

•  The  front  end  inductance  of  a  PRS  load  per  se  is  not  a  problem  for  energy 
coupling,  provided  this  inductance  remains  in  the  range  <20nH. 

•  The  gas  puff  PRS  load  on  Saturn  should  draw  at  least  8— *^10  MA  in  the  present 
configuration. 

The  observed  behavior^*  on  Saturn  -  10  MA  into  gas  puffs  and  12.5  MA  into 
short  circuit  loads  -  is  quite  nicely  in  line  with  this  prediction. 

Furthermore,  all  the  experience  to  date  with  the  more  complex  X-pinch  load 
design  indicates  that  (i.)  it  represents  a  very  credible  PRS  concept  in  terms  of  energy 
coupling  and  (ii.)  it  should  be  quite  compatible  with  the  presently  contemplated 
lOx  machine  designs.  The  model  we  have  discussed  here  will  be  refined  further  by 
2D  MHD  studies  to  include  the  role  of  runaway  electron  populations  in  the  energy 
budget  and  to  optimize  the  load  dynamics. 

More  general  problems,  involving  loss  mechanisms  in  MITL  segments,  can  be 
readily  addressed  by  the  method  and  some  initial  work  is  already  underway. 
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